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Abstract. In this paper we address the classical question going back to S. Bochner and 
H. L. Krall (see original [B], [Kr] and [EKLW] for the modern status quo) to describe all systems 
{p„(a;)}^Q of orthogonal polynomials (OPS) which are the eigenfunctions of some finite order 
differential operator, i.e. satisfy the equation 

N 

(1) ^afe(a;)y('=)(x) = A„y(a;). 

k=l 

Such systems of orthogonal polynomials are called Bochner-Krall OPS (or BKS for short) and 
their spectral differential operators are accordingly called Bochner-Krall operators (or BK- 
operators for short). We say that a BKS has compact type if it is orthogonal with respect 
to a compactly supported positive measure on M. It is well-known that the order A'^ of any 
BK-operator should be even and every coefficient {x) must be a polynomial of degree at most 
fc, see e.g. [EKLW]. Below we show that the leading coefficient of a compact type BK-operator 
is of the form {{x — a){x — h))^/"^. This settles the special case of the general conjecture 7.3. of 
[EKLW] describing the leading terms of all BK-operators. New results on the asymptotic distri- 
bution of zeros of polynomial eigenfunctions for a spectral problem (1) are the main ingredient 
in the proofs, comp. [BR] and [MS]. 



§1. Introduction 

Let PiK and Pc denote the spaces of all real and complex polynomials, respectively, in a 
variable x. By a real (or complex) polynomial system we will mean a sequence {Pn}'^=Q of 
polynomials in (or Pc) such that degp^ = n. By an orthogonal polynomial system (OPS) 
one understands a real polynomial system {pn} such that {pmPm) is nonzero precisely when 
n = m, where ( , ) is some reasonable inner product on the linear space Pr. If an orthogonal 
polynomial system for a given inner product exists, the pn are unique, up to multiplication 
by scalars. 

Orthogonal polynomial systems have been studied in various degrees of generality. Classi- 
cally, one has considered inner products of the form 

(p, q) = a{p ■ q) 
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where cr is a moment functional, that is a hnear functional on Pr. It is known that all moment 
functionals can be represented by an integral 



where /i is a (possibly signed) Borel measure on the real line. The most complete theories 
have been obtained in the case where /j, is positive and has compact support. In what follows 
we will mainly be concerned with orthogonal polynomial systems of this particular kind, 
which we will call compact type orthogonal systems . Recently, there has been interest in 
more general inner products called Sobolev, which are of the form 



where the cr/. are moment functionals. For the basics of the classical theory of orthogonal 
polynomials see e.g. [Sz] and [Ch]. 

Consider now a differential operator 



where the coefficients ak{x) are polynomials in Pc- We are interested in eigenpolynomials of 
this operator, that is polynomials p G Pc satisfying Dp — Xp for some constant A. Already 
S. Bochner observed that the operator D has infinitely many linearly independent eigenpoly- 
nomials if and only if dega^ < k, with equality for at least one k. In this case there is 
precisely one monic degree n eigenpolynomial p„ for all sufficiently large n. For generic aj-, 
the same is true for every n > 0. If a complex polynomial system consists of eigenpolynomials 
for an operator of the form (2), we will call it a system of eigenpolynomials. 

A Bochner- Krall system (BKS for short) is defined to be a real polynomial system which is 
both orthogonal (with respect to some inner product ( , )) and a system of eigenpolynomials 
(for some differential operator Xi). In this case ( , ) is called a Bochner- Krall inner product, 
and X) is called a Bochner-Krall operator. If a BKS is a compact type orthogonal system, we 
will call it for short a compact type BKS. The results we report in this note are valid for all 
compact type BKS. 

It is an open problem to classify all Bochner-Krall systems. A complete classification is 
only known for Bochner-Krall operators d with < 4. The corresponding BKS arc various 
classical systems such as the Jacobi-type, the Laguerre-type, the Legendre type, the Bessel 
and the Hermite polynomials, see [EKLW], Th. 3.1. In general, it is not even known which 
differential operators are Bochner-Krall operators for some BKS. In [EKLW] it is conjectured 
that the leading coefficient ajv of a Bochner-Krall operator is a power of a polynomial of 
degree at most 2. Our main result is an affirmative answer to this conjecture for compact 
type BKS. 




M 



(2) 




Our results are obtained by studying the asymptotic distribution of zeros of a polynomial 
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let CKi, . . . ,an denote the (real or complex) zeros of pn- Let Vn = -^Yl7=i^(^ ~ '^*) 
the probability measure in the complex plane with point masses at these zeros. We call 
the measures the root measures of the polynomial system {pn}- If the sequence of root 
measures converges weakly to a measure u when n — > oo, we say that u is the asymptotic 
distribution of zeros of the polynomial system. 

The following results, which characterize the asymptotic distribution of zeros for compact 
type OPS and for systems of eigenpolynomials respectively, are crucial to our treatment. 

Suppose that a polynomial system {pn} is orthogonal with respect to a positive measure 
/i, and that the convex hull of supp is a compact interval [a, b]. It is well known (see [Ch]) 
that the zeros of every pn are contained in the interval [a,b]. The following is a more precise 
result on the distribution of zeros. 

Theorem A, see [Ne], Th. 3, p. 50. Let the polynomial system {pn} be orthogonal 
with respect to a compactly supported, positive measure on M, and let the convex hull of 
supp/U be [a, b]. Then the asymptotic distribution of zeros of {pn} is an absolutely continuous 
measure u which depends only on [a,b]. The support of u is precisely [a, b] and its density in 
this interval is given by 

^ TTy/ (b — x){x — a) 

Next we describe the asymptotic distribution of zeros for a system of eigenpolynomials. 

Theorem B, see [BR] Th. 2 and 4. Let {pn} be a system of eigenpolynomials for an 
operator with ayv monic of degree A^. Then the asymptotic distribution of zeros of {pn} is 
a probability measure v with the following properties: 

a) V has compact support; 

b) its Cauchy transform C{x) = J satisfies the equation C{x)^ = l/aN(x) for almost 
all X E C 

These properties determine u uniquely. 

Note that the limiting measure v is independent of all terms in (2) except the leading term 

To derive from these two results a statement about compact type BKS, we will need the 
following. 

Proposition 1. Let {pn} be a system of eigenpolynomials for a differential operator 0, 
and assume that all the zeros of pn are real. Then there exists a compact set containing all 
the zeros of every pn if and only if degaAr = N. 

Now it is easy to derive the following 

Main Theorem. Let {pn} be a compact type BKS, orthogonal with respect to a measure 
/i and with differential operator D. If the convex hull of supp fj, is the interval [a, b], then N 
is even and aN{x) is a constant multiple of {{x — a){x — 6))^/^. 

References and acknowledgements. There exists a really vast literature devoted 
to the classification problem for OPS. Classification of BKS has also attracted considerable 
attention, see e.g. [EKLW] with its 100 references and [KLY1],[KLY2], [KYY] and refer- 
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active area. We are grateful to Dr. M. Shapiro and Professor H. Shapiro for a number of 
discussions on the topic. The third author wants to acknowledge the hospitality of Mathe- 
matisches Forschungsinstitut Oberwolfach in September 2001 in whose peaceful and serene 
atmosphere he found some information on BKS and realized that the results of [MS] and 
[BR] are applicable to the BKS-classification problem. 

§2. Proofs 

We need to prove Proposition 1 and the main theorem (as its corollary). Since the proof 
of Theorem B, in the situation where we will need it, follows easily along the same lines, we 
will include such a proof for the convenience of the reader. 

Consider a polynomial system {pn} with the associated root measures u^. Assume that 
the supports of the measures are all contained in the same compact set, and that ^ i' 
in the weak topology. Let Cn{x) be the Cauchy transform of Un and note that 

Cnix) - ^^^^^ 



npn{x)' 

If C{x) denotes the Cauchy transform of v, it follows that 



npn{x) 

for almost every a; e C. 



C{x) 



Suppose now that {pn\ is a system of eigenpolynomials for an operator 13 and that supp 
are all contained in the same compact subset of the real line. Then there exists at least a 
subsequence of the converging weakly to some measure u. Moreover, if we let denote 
the root measure of the jth derivative of then it follows from RoUe's theorem that (a 
subsequence of) converges weakly to v for every j > 0. In particular. 



(n - j)pn{x) 

for almost every x G C, where C{x) is the Cauchy transform of u. If we divide both sides of 
the differential equation dpn = ^nPn hy n{n — 1) ... {n — N + l)pn we obtain 



/=o {n - 3)p^d\x) n-N + 1 _ j)p\i\x) ' ' ' ~ n(n - 1) . . . (n - AT + 1) ' 



When n ^ oo all the terms in the left-hand side but the first tend to zero, and so 



aN{x)C{x)^ = lim 



n-»oo n{n — 1) ... (n — A" + 1) 

for almost all x. But it can be seen (sec [BR]) that = J2k=i (^k'n{n — 1) . . . (n — /c + 1) 
where is the coefficient at x^ in ak{x). In particular, if dega^v < A^ i.e. cat = 0, then 

\ I /n ( /n 1 ^ ( /n A/^ _l_ 1 ^ i. O tttVi on n i. r\/~\ o m tJ- ^r\ \ 1 rwiTd 
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X. This implies that u = 0, a contradiction. This argument proves one of the imphcations 
in Proposition 1. Moreover, if aAr(x) is monic of degree N, then cat = 1 so it foUows that 
C{x)^ = For the other imphcation we refer to [BR], Lemma 9. 

□ 

Suppose now that {pri{x)} is a compact type BKS as in the Main Theorem. By the remark 
preceeding Theorem A, the zeros of every Pn{x) are contained in the interval [a, b]. It follows 
from Proposition 1 that degaAr(x) = A^, and we might as well assume that aAr(x) is monic. 
Hence the Cauchy transform C{x) of the limit u = liuin^oo^n satisfies C{x)^ = l/aAr(x). 
On the other hand, a direct computation of the Cauchy transform, using the expression for 
u in Theorem A, gives C(x)^ = l/{x — a){x — b). Comparing these results yields aN{x) — 
((x-a)(x-6))^/2. 
□ 



§3. Final remarks 

Problem 1. The major problem in the context of this paper is whether every compact 
type BKS is a Jacobi-type OPS, compare [EKLW], Conjecture 7.3. For the constant leading 
coefficient the analogous fact was proved in [KYY]. 

Problem 2. Is there an analog of Theorem A on the asymptotic zero distribution for a 
signed measure n with compact support? What is the situation for a probability measure 
with a noncompact support as well as for Sobolev orthogonal polynomial systems. (There 
exists a literature on this topic.) 

Problem 3. Generalize the results of [BR] to operators with dega^ < N. Preliminary 
computer experiments show the existence of limiting measures for all operators (2). 
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